Kronos: a workflow assembler for genome analytics and informatics.

BackgroundThe field of next-generation sequencing informatics has matured to a point where algorithmic advances in sequence alignment and individual feature detection methods have stabilized. Practical and robust implementation of complex analytical workflows (where such tools are structured into "best practices" for automated analysis of next-generation sequencing datasets) still requires significant programming investment and expertise.ResultsWe present Kronos, a software platform for facilitating the development and execution of modular, auditable, and distributable bioinformatics workflows. Kronos obviates the need for explicit coding of workflows by compiling a text configuration file into executable Python applications. Making analysis modules would still require programming. The framework of each workflow includes a run manager to execute the encoded workflows locally (or on a cluster or cloud), parallelize tasks, and log all runtime events. The resulting workflows are highly modular and configurable by construction, facilitating flexible and extensible meta-applications that can be modified easily through configuration file editing. The workflows are fully encoded for ease of distribution and can be instantiated on external systems, a step toward reproducible research and comparative analyses. We introduce a framework for building Kronos components that function as shareable, modular nodes in Kronos workflows.ConclusionsThe Kronos platform provides a standard framework for developers to implement custom tools, reuse existing tools, and contribute to the community at large. Kronos is shipped with both Docker and Amazon Web Services Machine Images. It is free, open source, and available through the Python Package Index and at https://github.com/jtaghiyar/kronos

Feature-based classifiers for somatic mutation detection in tumour–normal paired sequencing data

Motivation: The study of cancer genomes now routinely involves using next-generation sequencing technology (NGS) to profile tumours for single nucleotide variant (SNV) somatic mutations. However, surprisingly few published bioinformatics methods exist for the specific purpose of identifying somatic mutations from NGS data and existing tools are often inaccurate, yielding intolerably high false prediction rates. As such, the computational problem of accurately inferring somatic mutations from paired tumour/normal NGS data remains an unsolved challenge

SNVMix: predicting single nucleotide variants from next-generation sequencing of tumors

Motivation: Next-generation sequencing (NGS) has enabled whole genome and transcriptome single nucleotide variant (SNV) discovery in cancer. NGS produces millions of short sequence reads that, once aligned to a reference genome sequence, can be interpreted for the presence of SNVs. Although tools exist for SNV discovery from NGS data, none are specifically suited to work with data from tumors, where altered ploidy and tumor cellularity impact the statistical expectations of SNV discovery

The driver landscape of sporadic chordoma.

Chordoma is a malignant, often incurable bone tumour showing notochordal differentiation. Here, we defined the somatic driver landscape of 104 cases of sporadic chordoma. We reveal somatic duplications of the notochordal transcription factor brachyury (T) in up to 27% of cases. These variants recapitulate the rearrangement architecture of the pathogenic germline duplications of T that underlie familial chordoma. In addition, we find potentially clinically actionable PI3K signalling mutations in 16% of cases. Intriguingly, one of the most frequently altered genes, mutated exclusively by inactivating mutation, was LYST (10%), which may represent a novel cancer gene in chordoma.Chordoma is a rare often incurable malignant bone tumour. Here, the authors investigate driver mutations of sporadic chordoma in 104 cases, revealing duplications in notochordal transcription factor brachyury (T), PI3K signalling mutations, and mutations in LYST, a potential novel cancer gene in chordoma

Multifocal clonal evolution characterized using circulating tumour DNA in a case of metastatic breast cancer.

Circulating tumour DNA analysis can be used to track tumour burden and analyse cancer genomes non-invasively but the extent to which it represents metastatic heterogeneity is unknown. Here we follow a patient with metastatic ER-positive and HER2-positive breast cancer receiving two lines of targeted therapy over 3 years. We characterize genomic architecture and infer clonal evolution in eight tumour biopsies and nine plasma samples collected over 1,193 days of clinical follow-up using exome and targeted amplicon sequencing. Mutation levels in the plasma samples reflect the clonal hierarchy inferred from sequencing of tumour biopsies. Serial changes in circulating levels of sub-clonal private mutations correlate with different treatment responses between metastatic sites. This comparison of biopsy and plasma samples in a single patient with metastatic breast cancer shows that circulating tumour DNA can allow real-time sampling of multifocal clonal evolution.We thank the Human Research Tissue Bank at Addenbrooke’s Hospital which is supported by the NIHR Cambridge Biomedical Research Centre. We acknowledge the support of Cancer Research UK, the University of Cambridge, National Institute for Health Research Cambridge Biomedical Research Centre and Cambridge Experimental Cancer Medicine Centre. Dr. Dawson was supported by an Australian National Breast Cancer Foundation and Victorian Cancer Agency Early Career Fellowship. Dr. Murtaza was supported by Science Foundation Arizona’s Bisgrove Scholars Early Tenure Track award.This is the final version of the article. It first appeared from Nature Publishing Group via http://dx.doi.org/10.1038/ncomms976

Model-based clustering of array CGH data

Motivation: Analysis of array comparative genomic hybridization (aCGH) data for recurrent DNA copy number alterations from a cohort of patients can yield distinct sets of molecular signatures or profiles. This can be due to the presence of heterogeneous cancer subtypes within a supposedly homogeneous population

Regulation of pH by Carbonic Anhydrase 9 Mediates Survival of Pancreatic Cancer Cells With Activated KRAS in Response to Hypoxia.

Background & Aims Most pancreatic ductal adenocarcinomas (PDACs) express an activated form of KRAS, become hypoxic and dysplastic, and are refractory to chemo and radiation therapies. To survive in the hypoxic environment, PDAC cells upregulate enzymes and transporters involved in pH regulation, including the extracellular facing carbonic anhydrase 9 (CA9). We evaluated the effect of blocking CA9, in combination with administration of gemcitabine, in mouse models of pancreatic cancer. Methods We knocked down expression of KRAS in human (PK-8 and PK-1) PDAC cells with small hairpin RNAs. Human and mouse (KrasG12D/Pdx1-Cre/Tp53/RosaYFP) PDAC cells were incubated with inhibitors of MEK (trametinib) or extracellular signal-regulated kinase (ERK), and some cells were cultured under hypoxic conditions. We measured levels and stability of the hypoxia-inducible factor 1 subunit alpha (HIF1A), endothelial PAS domain 1 protein (EPAS1, also called HIF2A), CA9, solute carrier family 16 member 4 (SLC16A4, also called MCT4), and SLC2A1 (also called GLUT1) by immunoblot analyses. We analyzed intracellular pH (pHi) and extracellular metabolic flux. We knocked down expression of CA9 in PDAC cells, or inhibited CA9 with SLC-0111, incubated them with gemcitabine, and assessed pHi, metabolic flux, and cytotoxicity under normoxic and hypoxic conditions. Cells were also injected into either immune-compromised or immune-competent mice and growth of xenograft tumors was assessed. Tumor fragments derived from patients with PDAC were surgically ligated to the pancreas of mice and the growth of tumors was assessed. We performed tissue microarray analyses of 205 human PDAC samples to measure levels of CA9 and associated expression of genes that regulate hypoxia with outcomes of patients using the Cancer Genome Atlas database. Results Under hypoxic conditions, PDAC cells had increased levels of HIF1A and HIF2A, upregulated expression of CA9, and activated glycolysis. Knockdown of KRAS in PDAC cells, or incubation with trametinib, reduced the posttranscriptional stabilization of HIF1A and HIF2A, upregulation of CA9, pHi, and glycolysis in response to hypoxia. CA9 was expressed by 66% of PDAC samples analyzed; high expression of genes associated with metabolic adaptation to hypoxia, including CA9, correlated with significantly reduced survival times of patients. Knockdown or pharmacologic inhibition of CA9 in PDAC cells significantly reduced pHi in cells under hypoxic conditions, decreased gemcitabine-induced glycolysis, and increased their sensitivity to gemcitabine. PDAC cells with knockdown of CA9 formed smaller xenograft tumors in mice, and injection of gemcitabine inhibited tumor growth and significantly increased survival times of mice. In mice with xenograft tumors grown from human PDAC cells, oral administration of SLC-0111 and injection of gemcitabine increased intratumor acidosis and increased cell death. These tumors, and tumors grown from PDAC patient-derived tumor fragments, grew more slowly than xenograft tumors in mice given control agents, resulting in longer survival times. In KrasG12D/Pdx1-Cre/Tp53/RosaYFP genetically modified mice, oral administration of SLC-0111 and injection of gemcitabine reduced numbers of B cells in tumors. Conclusions In response to hypoxia, PDAC cells that express activated KRAS increase expression of CA9, via stabilization of HIF1A and HIF2A, to regulate pH and glycolysis. Disruption of this pathway slows growth of PDAC xenograft tumors in mice and might be developed for treatment of pancreatic cancer

A tumor DNA complex aberration index is an independent predictor of survival in breast and ovarian cancer

Complex focal chromosomal rearrangements in cancer genomes, also called "firestorms", can be scored from DNA copy number data. The complex arm-wise aberration index (CAAI) is a score that captures DNA copy number alterations that appear as focal complex events in tumors, and has potential prognostic value in breast cancer. This study aimed to validate this DNA-based prognostic index in breast cancer and test for the first time its potential prognostic value in ovarian cancer. Copy number alteration (CNA) data from 1950 breast carcinomas (METABRIC cohort) and 508 high-grade serous ovarian carcinomas (TCGA dataset) were analyzed. Cases were classified as CAAI positive if at least one complex focal event was scored. Complex alterations were frequently localized on chromosome 8p (n = 159), 17q (n = 176) and 11q (n = 251). CAAI events on 11q were most frequent in estrogen receptor positive (ER+) cases and on 17q in estrogen receptor negative (ER) cases. We found only a modest correlation between CAAI and the overall rate of genomic instability (GII) and number of breakpoints (r = 0.27 and r = 0.42, p <0.001). Breast cancer specific survival (BCSS), overall survival (OS) and ovarian cancer progression free survival (PUS) were used as clinical end points in Cox proportional hazard model survival analyses. CAAI positive breast cancers (43%) had higher mortality: hazard ratio (HR) of 1.94 (95%CI, 1.62-2.32) for BCSS, and of 1.49 (95%CI, 1.30-1.71) for OS. Representations of the 70-gene and the 21-gene predictors were compared with CAAI in multivariable models and CAAI was independently significant with a Cox adjusted HR of 1.56 (95%CI, 1.23-1.99) for ER+ and 1.55 (95%CI, 1.11-2.18) for ER disease. None of the expression-based predictors were prognostic in the ER subset. We found that a model including CAM and the two expression-based prognostic signatures outperformed a model including the 21-gene and 70-gene signatures but excluding CAAL Inclusion of CAAI in the clinical prognostication tool PREDICT significantly improved its performance. CAAI positive ovarian cancers (52%) also had worse prognosis: HRs of 1.3 (95%CI, 1.1-1.7) for PFS and 1.3 (95%CI, 1.1-1.6) for OS. This study validates CAM as an independent predictor of survival in both ER+ and ER breast cancer and reveals a significant prognostic value for CAAI in high-grade serous ovarian cancer. (C) 2014 The Authors. Published by Elsevier B.V. on behalf of Federation of European Biochemical Societies. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).Publisher PDFPeer reviewe